Liquid Crystal Analogue of Abrikosov Vortex Flow in Superconductors 
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' 

We extend the correspondence between the Renn-Lubensky Twist-Grain-Boundary-A phase in 
chiral liquid crystals and the Abrikosov mixed state in superconductors to dynamical aspects. We 
find that for a TGB sample with free boundaries, an external electric field applied along the helical 
axis induces a uniform translational motion of the grain boundary system - an analogue of the well- 
known mixed state flux flow. Likewise, an analogue of the mixed state Nernst effect is found. In 
much the same way in which the flux flow carries intercore electric fields generating Joule heat in an 
otherwise dissipation-free system, the grain boundary flow carries along polarized charges, resulting 
in a finite electric conductivity in a ferroelectric. 

PACS numbers: 61.30.Cz, 61.30.Jf 



The formal analogy of the smectic free energy within a 
Ginzburg-Landau approach, with that for the well known 
case of a superconductor has had a long history in the 
study of liquid crystals. It was first pointed out by 
de Gennes that the essential mean field features of the 
phase transition between the Nematic (N) and Smectic- 
A (SmA) phases can be derived by simply resorting to the 
correspondence with the Normal-Superconductor transi- 
tion. Q Corrections to this mean-field picture was stud- 
ied by means of a large-N expansion of the GL energy 
functional, in which a runaway RG trajectory was inter- 
preted as an indication of a fluctuation-driven first order 
transition. Q 

An important application of this analogy was carried 
out by Renn and Lubensky, who, based on the mean- 
field phase diagram of type-II superconductors in an ap- 
plied magnetic field, predicted an intermediate phase 
lying between the cholesteric (Ch) and SmA phases. 
H This new phase, termed the Twist-Grain-Boundary- 
A(TGBA) phase, is the counterpart of the Abrikosov vor- 
tex lattice phase (also known as the mixed phase) , and is 
composed of a regular network of screw dislocations. For 
the detailed structure of this lattice, readers are reffered 
to ref . H . The TGB A phase was also independently dis- 
covered by exerimenters and immediately became a 
subject of intense research. 

So far, theoretical investigations on the TGBA based 
on the superconductor analogy have been restricted to 
mean-field, static aspects. (See however refs Con- 
sidering however the rich variety of mixed state physics 
that continues to be revealed || through recent (high 
Tc-motivated) studies of fluctuation effects and vortex 
dynamics, it is natural to ask what new feature arises 
when one seeks to extend the analogy in similar direc- 
tions. As in high Tc superconductors, it is expected that 
thermal fluctuations have profound effects on the TGBA 
lattice structure; a simple estimate via the Ginzburg cri- 
teria strongly suggests this vulnerability, and there is an 
additional intrinsic instability towards disorder which is a 



consequence of the very symmetry of the TGB. How- 
ever, a detailed study of the full effect of thermal fluc- 
tuation within the GL formalism requires considerable 
computing efforts, and we choose to set a less ambitious 
goal. In this paper, we study the fundamental properties 
of dislocation motion in the TGB system. As a result, we 
obtain several new entries to the list of superconductor- 
liquid crystal analogue, as sumarised in Table 1. A major 
feature is the precise liquid crystal analogue of the flux 
flow motion of the Abrikisov vortex lattice. Although 
the long wavelength hydrodynamics of the TGBA phase 
has been previously worked out, [fiof this paper is, to 
our knowledge, the first one which has considered the 
grain boundary itself and its constituent screw disloca- 
tion lines as mobile entities. An important implication 
of this motion is that the dislocation cores carry elec- 
tric charges along with the flow through the ferroelectric 
media. Thus the system is predicted to be conductive. 

To generate the motion of the dislocations, we will take 
advantage of the chiral nature of system. Due to the 
lack of chiral symmetry, polar vector fields can couple 
linearly to axial vector fields. If we consider the torque 
of the smectic layer normal as the latter and couple it to 
an external polar vector field, this will effectively exert 
a force on the grain boundary. As the external field, we 
consider (1) an electric field applied along the helical axis 
of TGBA phase and (2) temperature gradient imposed 
along the same direction. 

In the SmA phase of such materials, it is well known 
that when an electric field is applied in-plane to the smec- 
tic layers, a finite tilting of the molecular axis from the 
smectic layer normal is favored. This phenomenon is 
called the electroclinic effect, jll] We shall show here that 
the stress induced by the electroclinic effect produces the 
driving force acting on screw dislocations. Our starting 
point is the following free energy, which is a combination 
of the covariant Chen-Lubensky(CL) model and the 
energy related to the electroclinic effect. 
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where, n is the Frank director, ^(x) is the smectic order 



parameter, ip(x) 
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with go = 2ir/d, d 



the separation between smectic layers, u the smectic layer 
displacement, n is the normal vector of smectic layer, -q 
is the coupling constant between the polarization P and 
n x (Vit — (5n), and \ is the electric susceptibility. The 
terms ?7(no x (Vit — <5n)) • P and Kik^w • V x n reflect 
the mentioned symmetry of the system. Differentiating 
F with respect to <5n, gives the stress field J = — = 

B'(y - Sn) with v = Vu, and B' = 2Cq%\ $ f - rfx- 
J is analogous to the supercurrent 3 S which is the free 
energy differentiated with respect to the vector potential 
A. x-axis, ^From the solutions of V ■ J = and 4^ = 
0, we see that v — 5n v — <5n e is divided into two parts. 
One is the "vortex" part v — Sn v , 12 with nonvanishing 
curl, and the other is the uniform part, <5n e = T,x ^, xE ^ , 
induced by the electric field E. Applying the virtual work 
principle, we find a Peach-Koehler tyoe force 



f/ = J e x di 



(2) 



acting on each unit length of the dislocation. Here di is 
the Burgers vector. 

Now we compare this with superconductors. B'(y — 
6n v ) corresponds to the vortex part of the supercurrent, 
and — B'6n e corresponds to the transport current. In 
the superconductor case, the corresponding force fx, = 
Jt x <£>oZ is just the Lorentz force acting on the vortex 
line per unit length. $o is the flux quantum. Compar- 
ing f; with fx, it is clear that the force acting on the 
single screw dislocation corresponds to Lorentz force in 
superconductors. 

Next we consider the case where the Peach-Koehler 
force acts on the grain boundaries which is just an array 
of screw dislocations. For the grain boundary located at 
x = 0, the unit normals rip and of the SmA slab at 
x < and x > , respectively are forced to be rotated 
with respect to each other around the x-axis by an angle 
of 2na. A uniform electric field E = E~k applied along the 
x-axis will generate the stress "currents" on the left and 
right ends of each SmA slab given by Jg = — ?7X n o x E, 
and Jf- = —rix^-Q' x E, respectively (See Figl). The 
extension of eq.(gj) to this configuration is straightforward 
and the force itgb acting on the grain boundary per unit 
area is 



Itgb 



r/xEd 



cos(7ra)x ~ —2i]xETra£, 



(3) 



where the final expression is valid for small a. This sim- 
ple superposition is invalidated near k = k c i where dis- 
locations cores begin to overlap. 



Next we consider the case where a temperature gra- 
dient is imposed along the pitch. In analogy with the 
Nernst effect in a superconducting mixed phase, it is 
anticipated that screw dislocations in the TGBA phase 
j^a+rc influenced by a thermal force generated in the di- 
rection of the temperature gradient. We shall show below 
that this is indeed the case. The physics of what is shown 
below may be summarized as follows. The temperature 
gradient produces the gradient of smectic density | i/> |, 
where n s = \ ip \ . This in turn disturbs the equilibrium 
of screw dislocations, giving rise to an effective driving 
force acting on the system. 

Following the approach of Koklov and co-workers |D| 
who evaluated the Nernst coefficient for the superconduc- 
tor case solely within the Ginzburg-Landau formalism, 
we can derive Q an expression for the total force acting 
on the ensemble of screw dislocations in the case where 
the temperature gradient produces a steady nonequilib- 
rium distribution of screw dislocations. For the condition 
(i.e. £ s 2 /hid <C 1 with £ s the smectic correlation length), 
the effective force per unit area on the grain boundary is 



f T = A 
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B = 2Cq$ \ ip | , arid the penetration length A2 = 
\J K2/B. It should be noted that the drift force f is 
proportional to the gradient of smectic density V(An s ), 
The first term of the right hand side of eq. (|^) has simple 
form of a single dislocation contribution multiplied by a 
factor of I /Id, the density of screw dislocations in a layer. 
The second term reflects the global twisted structure of 
the TGBA system. 

Now that we have obtained the forces induced by the 
coupling to external fields, we now turn to the actual 
motion of the dislocations. 

We start with the "model A" type equation for the 
order parameter (Til : 



5F 

= + C (t) 



(5) 



((t) is a thermal noise obeying the fluctuation-dissipation 



is 



relation, < C(*)C(f) >= f^S(t - if). ^ = ^ + r$ 
generically a complex kinetic coefficient. Once the defect 
moves, we will have to consider the effect of a hydrody- 
namic lift force (Magnus force). The existence of such 
a force and its relevance to the Hall sign anomaly is an 
issue of present debate in superconductors. Jll|[ll| It can 
be shown that a imaginary part of the relaxation time 
7^ in cq. (refptd) can give rise to such a force. How- 
ever, renormarization-group calculation shows (7^/7^,) 
term to be an irrelevant perturbation in the vicinity of 
the transition point, i. e. near k — k c2 . Furthermore, 
it was concluded in Ref. JTsJ that the force is absent in 
the single screw dislocation case. Thus we consider it 
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plausible to neglect such forces, and put 7^/7^ = 0. In 
this purely dissipative case, the motion of screw disloca- 
tion takes the form 7;v; = f;, where 7; is the viscosity 
coefficient for a screw dislocation and v; is its velocity. 
7; is determined using the Rayleigh functional for eq.(g) 
R[d t ip,d t tp*] = J d 3 x^ I d t tp | 2 , and the balance if force 
= f). as 7; = 7^7r ln(A/£ s ). Hence the weak chi- 
rality regime, eq.(||) reduces to the stochastic equation 
for a single defect line, 7;^- = — jjjf + Gj(t). The nois 
Gi(t) satishes < Gi(t)Gi(t') >= ^^8(t-t'). Later we 

will show that the velocity vi = T, ^' E is obtained from 
the experimental parameters 77, Xi d. 

At ko ~ k C 2 {hid ~ £ 2 ) f ne screw dislocations cannot 
be treated as independent mobile entities. The static so- 
lutions of the linearized GL equation for if) in the super- 
conductor case near H C 2 can be approximated by eigen- 
functions of the lowest Landau level. An analogous pro- 
cedure can be taken for the TGBA case when qo/lo ^S> 1. 
P| Here we add the time dependence into the solution. 



*(*) 



exp 



C s exp (iq s ± (x s + U s (t)) ■ x 
(k x - ko( Xj: + U s (t))) 2 

e 



(6) 



Here, x s is the position of s-th SmA slab, U s (t) 
denotes the displacement of the zero of the s-th 
eigenfuntion in the lowest Landau level. qj_(x s + 
U s (t)) = q o (0,sm(k o {x s + U s {t)),cos{k (x s + U s (t))), 
£ 2 = l/(qoko). The dissipation function R per unit vol- 



ume in this case R 



hqa 



. As we will discuss next, 



the velocity vi (now defined as the average velocites of the 
zeroes) is decided by combining R and the electric energy 
dissipation. 

First we point out that a d.c. polarization current is 
directly observable as a direct consequence of the mo- 
tion of screw dislocations. This is in sharp contrast to 
dislocation-free feroelectrics, where only a.c. currents are 
allowed to flow. The origin of the charge carried by the 
screw dislocation stems from the chiral nature of the sys- 
tem. As can be seen from eq.(l), the deviation of the 
director vector from the layer normal gives rise to polar- 
ization. Around a screw dislocation this is obtained as 



■/Co ( ^ X a^' ' ) e r + X e oE per unit length, where 
/Cq is the modified Bessel function, e r = (cos 4>, sin<^>, 0), 



rjd 
' 2ttX 



and d> = tan 



. Hence at the surface of the core 



of a screw dislocation we have a local charge accumula- 
tion of Q = J e?SV ■ P = f)xd- Thus near k c \ we are 
able to relate the velocity to the polarization current by 
summing up the contribution dtQ and get 



rjxvikoxx. 



(7) 



Here, ko is the spatialy averaged twist of the director, for 
which the relation fco = j4~ f° r small a is used. Near k C 2 



(VWb ~ £s ) using the solution of eq. (||) leads to the 
same form as in eq. (Q). As is the case of superconductors 
fco is expressed in the form of the expansion of 



fcc2 — fco 
fc 



near 



k C 2- The parameter (3 [|3j (corresponding to Abrikosov's 
beta in superconductors) is related to the configuration of 
screw dislocation lattice. Differentiating the free energy 
with respect to ^fco, we obtain fco = fco — fc ^ 2 7^° with 



«2 = J ^y 2 - . Demanding the dissipation due to the 
polarization current to coincide with R which is induced 
by the defects dynamics, we are able to express vi as 
in terms of the material parameters fco, 



7v k l d ^/Tr 



VI = 

77, Xj ancl viscosity7^,. 

In summary we have investigated the basic electric 
properties of the dislocation motion in the TGBA phase, 
with a special emphasis on the analogy with vortex dy- 
namics in superconductors. In actual experimental sit- 
uations, rubbing process to fix the direction of directors 
of molecules should not be made in order to make the 
predicted effects observable. The present work should 
be considered a starting point for further studies of the 
dynamics in the TGB system, in which the interplay 
between thermal fluctuation and defect dynamics is ex- 
pected to give rise to various interesting physics. 

It is pleasure to acknowledge many helpful discussions 
with Professor T. Tokihiro at the early stage of this work. 
A. T. acknowledge a fellowship from the Japan Atomic 
Energy Research Institute. T. O. was supported in part 
by the Japan Society for the Promotion of Science. 



[1] P. G. de Gennes, Solid State Commun, 14, 997 (1973). 
[2] B.I.Halperin, T. C. Lubensky, and S. K. Ma, Phys. Rev. 

Lett., 32, 292 (1974). 
[3] S. R. Renn and T. C. Lubensky, Phys. Rev. A, 38, 2132 

(1988). 

[4] J. Goodby et. al, Nature (London) 337, 449 (1988). 
[5] T. Ota, A. Tanaka and R. Hayakawa, unpublished. 
[6] K. S. Bedell, M. Inui, D. Meltzer, J. R. Schrieffer, and 
S. Doniach, Phenomenology and Applications of High- 
Temperature Superconductors (Addison Wesley, 1992). 
[7] J. Toner, Phys. Rev. B, 43, 8289 (1991). 
[8] T. C. Lubensky, T. Tokihiro, and S. R. Renn, Phys. Rev. 

A, 43, 5449 (1991). 
[9] R. D. Kamien and T. C. Lubensky, J. Phys. I. France, 3, 
2131 (1993). 53, 650 (1996). 
[10] Y. Hatwalne, S. Ramasaway and J. Toner, Phys. Rev. 

Lett., 70, 2090 (1993). 
[11] P. S. Meyer, Mol. Cryst. Liquid Cryst., 40, 33 (1977). 
[12] A. R. Day, T. C. Lubensky, and A. J. McKane, Phys. 

Rev. A, 27, 1461 (1983). 
[13] A. Kuklov, A. Bulatov, and J. L. Birman, Phys. Rev. 
Lett., 72, 2274 (1994). 



3 



FIG. 1. The stress field near the grain boundary. J„ is the 
nonvanishing curl stress field around screw dislocations with 
its separeation Id- Je and is the electric field induced 
spatially uniform stress field of left SmA slab and right SmA 
slab, respectively. 

FIG. 2. The smectic layers near the grain boundary, d 
is layer separation of SmA slabs. The angle between electric 
field induced stress field in left SmA slab Jg along the Sm 
layer and that of right SmA slab is 1-kol. 
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superconductors 



- :supercurrent j 



liquid crystals 



-"^"■stress 

on 



Lorentz force 
vortex lattice flow 

Nernst effect 

flux flow resistivity 



Peach-Koehler force 

grain boundary flow 

due to the electric field 

grain boundary flow 

due to the temperature gradient 

conductivity due to grain boundary flow 



TABLE I. Dynamical aspects of the analogy between su- 
perconductors and liquid crystals. 
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